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Abstract 



The use of MS-like renormalization schemes in QCD requires an imple- 
mentation of nontrivial matching conditions across thresholds, a fact often 
overlooked in the literature. We shortly review the use of these matching con- 
ditions in QCD and check explicitly that the prediction for as{Mz), obtained 
by running the strong coupling constant from the M^- scale, does not substan- 
tially depend on the exact value of the matching point chosen in crossing the 
6-quark threshold when the appropriate matching conditions are taken into 
account. 
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During the last years a great effort has been done at LEP in order to measure the 
strong couphng constant as{Mz) at the Z mass scale [|l], H, ^. This measurement has 
been of crucial importance since it allowed, within the experimental errors, the running 
of the strong coupling constant to be checked from low energies to the electroweak scale. 
However, by going to higher orders in the renormalization group equations, some confusion 
has arisen in the literature on the different prescriptions one could use to cross thresholds 
in the evolution of the running coupling constant. The problem appears when working 
in MS*- like renormalization schemes: since these are mass-independent, the decoupling 
theorem of Appelquist-Carazzone is not fulfilled in "non-physical" quantities such as 
beta functions or coupling constants. Only in physical quantities particles with large 
masses do decouple. Logarithms of large masses induced by the renormalization group 
equations in the couplings are cancelled against other logarithms that appear in the 
calculation of physical observables. This is obviously an inconvenient, since a lot of effort 
must be invested in intermediate stages of a calculation to compute terms that will cancel 
in physical quantities. To remedy this problem the standard procedure has been the use 
of the effective field theory language P, 0, H]. For example, in QCD with a heavy quark 
and — 1 light quarks, one builds a theory with quarks and an effective field theory 
with — 1 quarks. Around the threshold of the heavy quark one requires agreement of 
the two theories. This gives a set of matching equations that relate the couplings of the 
theory with A^ quarks with the couplings of the theory with A^ — 1 quarks. This way, 
below the heavy quark threshold one can work with the effective theory, but using effective 
couplings. Then, by construction, decoupling is trivial. This procedure is equivalent to 
other renormalization schemes and allows us to correctly obtain the asymptotic value of 
the coupling constant. The price one has to pay is that coupling constants might not 
be continuous at thresholds. All this machinery is well established since the early 80's 
HI) B H' 1131 matching conditions were computed at the one-loop level |]^, || and 
at the two- loop level 10 1 for general gauge theories. It also seems to be well known 
for people working with GUTs where, in general, special attention has been paid to 
matching conditions at the different thresholds. However, the fact that one has to use 
appropriate matching conditions in passing thresholds has been frequently overlooked in 
the running of the QCD coupling constant by just taking a continuous coupling constant 
across thresholds. Then, the final results depend strongly on the exact scale one uses to 



connect the couplings 12, 13, H]. To solve this ambiguity some of the authors |11 



vary the matching scale between 0.75 and 2.5 times the mass of the heavy quark; others 
T2| , p!3| use directly Hth = '^niq, and yet others determine as(M^) with both fith = 



and fith = 2mq and take the average. Here we will show that when appropriate matching 
conditions are taken into account the final answer does not depend on the exact ^th used 
to connect the couplings. 

Although most of the points discussed in this paper are well known in some circles, 
given the confusion that exists in the literature and the importance of the subject we 
found it convenient to recall what the correct matching conditions are and to show that 
when they are consistently taken into account the dependence on the renormalization 
scale cancels (at least at the order the calculation is done). Consistency requires that 
if the evolution of the gauge coupling constant is done at n loops, matching conditions 
should be imposed using n — 1 loop formulae 0; then the residual dependence on the 
renormalization scale is of order n + 1. We will show this, explicitly, when running the 
QCD gauge constant from the r mass to the Z mass passing through the h threshold. 

The renormalization group equations in QCD for the strong gauge coupling constant 
and the quark masses are 
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where 

and /io is some reference point. 

The j3 coefficients governing the evolution of the gauge couphng constant are 

Pi = (4) 

^ 2 V 9 27 ^ 

with Np the number of quark flavours with mass lower than the renormalization scale /i. 
The first two coefficients are scheme-independent (in MS'-like schemes) but the higher- 



order coefficients depend on the renormalization conditions [|14|. We give P2 in the MS 
scheme. 

The quark mass anomalous dimensions are 

70 = 2 

71 = Np . (5) 

' 12 18 ^ ' 

Integration of eq. (|l|) can be performed by ffist inverting the series on the right-hand 
side of eq. (|l|) and then integrating on a and t. Finally, one can solve for a, at the required 
order, by using iterative methods. The result we obtain can be written in the following 
form 

asifi) = a^^\fi) + af\^i) + af\ii) + ■■■ (6) 

where a'f'\ af'^ represent the one-, two- and three- loop contributions respectively, 
and are given by 

= T-^^4H^ (7) 
1 + asifXoJPot 
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(/i) = ("i'^)' {bllogK{^^) (logir(/i) - 1) - [bl - h) (1 - K{f,))) , (9) 

where 

The running quark mass can also be obtained analytically at the one-loop level, from 
eq. (0). The result is 

m'ifi) = m2(^o)(l + a.(/io)/?ot)-'^"/^° • (H) 



For as{fJ.o)PQt « 1 the result is independent of Po and can be written as (we use 70 = 2) 



m\fi) = m\fio){l - 8a,(/io)t + ■■■)• (12) 

This expression can be used as long as fi is not very different from /iq; in particular we 
could use it to simplify the matching conditions. 



Conventionally ||T^, |I4[, higher-order RGEs are solved by doing a power series ex- 
pansion in 1/L with L = log(/i^/A^). This solution is given in terms of the QCD scale 
A, which is defined in such a way that it is renormalization- group-invariant but scheme- 
dependent and the so-called invariant mass rh. Passing of thresholds is implemented by 
requiring continuity of the couplings at threshold, which in turn requires defining different 
A's for different Np. For our purposes we prefer to use the solutions given above because 
they allow us to work more easily with scale-dependent matching conditions. 

In the MS scheme, or any of its simple modifications such as MS, the beta func- 
tion governing the running of the strong coupling constant is independent of quark 
masses. Then, contrary to what happens in momentum-subtraction schemes (MO), the 
Appelquist-Carazzone theorem ^ that states, when it can be applied, that the heavy par- 
ticles decouple at each order of perturbation theory is not realized in a trivial way. The 
decoupling of the heavy particles is fulfilled in physical quantities, but coupling constants 
and beta functions do not exhibit it. 

To obtain decoupling in MS schemes we need to build in the decoupling region, 
^ ^ M, an effective field theory that behaves as if only the light degrees of freedom 
were present. Matching conditions connect the parameters of the low-energy effective 
Lagrangian with the parameters of the full theory. This can be done by evaluating some 
Green functions in perturbation theory with both the full and the effective theories, then 
require they are the same, up to terms 0(1/M), for values of the renormalization scale 
just around the threshold. Then, the coupling constant of the effective theory can be 
expressed as a power series expansion in the coupling of the full theory with coefficients 
that depend on log(M//i). In order to obtain a good approximation using only the first 
few terms in the perturbative expansion, we have to evaluate matching conditions in a 
region where M/fi ~ 0(1). However, the results of these calculations should not depend 
on exactly which /j, is chosen. 

One-loop matching conditions have been obtained in |^ for a general gauge theory. 
To obtain matching conditions in QCD at the two-loop level several approaches have been 
pursued. Ovrut and Schnitzer computed the gluon self energies at the two-loop level 
with both the full and the effective theories and then required matching in the threshold 



region. We will follow a more direct approach, devised by Bernreuther and Wetzel ||10|| . 
Using the MO scheme as an intermediate stage, these authors were able to relate the 
MS coupling constant aj^{n), with Np quark flavours, with the gauge coupling constant 
a^j^{fi) of the effective fleld theory with Np — 1 quark flavours in which a heavy quark with 
TTijjg mass has been integrated out. This is because in momentum subtraction schemes 
the decoupling theorem is also realized in the coupling constants. The obtained relation 
has the following form: 



^MS = ^MS i'^ + Y.^JTS^kix)] (13) 

with 



\ k=l 



^ = ^M^W/^') (14) 



In order to calculate the coefficients Ck Bernreuther and Wetzel impose the RGEs, eq. (|T]) 
and eq. (@), on ajjg, and m-jjg, and obtain for the ffist two coefficients a set of 
coupled ffist-order linear differential equations depending only on the beta and gamma 
functions of the full and the effective theories. By solving them they found for a general 
SU (N) group the following result valid for the MS schemeQ 

;C2{R) + 777^ (^C2{G) - G2{R)] ^ Tr{I}\n=, 
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with C2{G) = N and C2{R) = — 27v — Casimir operator eigenvalues of the adjoint and 
fundamental representations, respectively, and D the space-time dimension. A technical 
point about the trace of the identity in the Dirac space, Tr{J}, should be discussed here. 
Strictly, in a general D-dimensional space, D even, the only irreducible representation of 

{^^,Y} = '^g"' (16) 

has dimension f{D) = 2-^/^. However, we can choose Tr{J} = f{D) = 4, or any other 
smooth function with /(4) = 4. Different choices of f{D) lead to different trivial modifi- 
cations of the MS renormalization scheme. However, as can be seen in eq. (|T5D , different 
choices give quite different matching conditions. Hence, in order to specify completely 
the renormalization scheme within the MS'-like schemes one should also specify which 
convention has been used for Tr{/}. Here we will use the usual convention among phe- 
nomenology papers, i.e. Tr{J} = 4. Then for QCD we have the following two-loop 
matching condition to connect the theory with — 1 quarks with the theory with N 
quarks at the g-quark threshold |T^ 

aN-iifJ'th) = OiNilJ'th) H 7: log 

OTT Hth 

iog!^V + ^iog!^ + i|. (17) 

97r lith ) 4 iith 

Here, ^th is the value at which we require matching. As commented, this equation is valid 
for arbitrary values of ixth as long as it is not far away from mq{mq). Should we use instead 
Tr{/} = 2-^/^, the logarithm in the second term would be changed to \og\/2mg/ nth 
changing completely the behaviour of the matching conditions. For example, it is clear 
from the above equation that one can always choose a fj,th in order to make the coupling 
continuous across thresholds. Using only the one-loop matching condition, i.e. only the 
ffist two terms in the right-hand side of eq. (|T^, and with Tr{/} = 4, we should require 
QN-ii^q) = CiNi^q)- Usiug the two-loop matching condition, the matching point is 
slightly different 0. However, if a scheme with Tr{/} = 2^/^ is used the matching point is 

^The solution of the two (for two loops) differential equations depends on two scheme-dependent 
arbitrary constants. To fix them one has to perform a complete calculation in the scheme one is interested 
in. 

^One can still impose aM-iirriq) — aN^mq) at the two-loop level, as Marciano does but this 
requires a slight modification of the MS scheme in order to absorb the non-logarithmic term in eq. (^) . 




A 



found around \/2mq. In what follows we will keep the couplings discontinuous and check 
the invariance of the final result with respect to the chosen matching scale ixth- 

Equation ( [T7| ) can be simplified by using eq. (|T^ to remove the dependence on the 
running mass and leave the result in terms of = mq{mq), the MS running mass 
evaluated at its own value. Since we are running from low energies to high energies it is 
also better to use the inverted equation. Thus, our matching condition at the threshold 
of the quark q will be: 



(^NifJ'th) = OiN-lifJ'th) log 



o log — 

37r nth 



+ ^?^i^fflog!^V-^log^-^l . (18) 

We start from a scale below the bottom-quark threshold, where we know the value of 
the strong coupling constant^ , e.g. Mr 

Mr = 1776.9 ± 0.7 MeV 

(19) 

asiMr) = 0.36 ± 0.03 , 

We use0 eq. ([18|) with fith = M^- and rriq = rric to obtain a^^^Mr) in terms of a^{Mr). 
Then we evolve a4(yu) until the Z boson mass scale by imposing matching conditions 
at an arbitrary intermediate scale ^th around m^. To run a4(/i) from M^. until ^th "we 
use eq. (|) with /io = Mr and four-quark beta functions. Then at ^th we impose the 
matching condition eq. ( pTS]) with ruq = mi, to obtain a^{^th) in terms of a^^fith)- Finally, 
to run a5(/x) from fith to Mz we use again eq. (^, but now with /xq = fith, and with 
five-quark beta functions. The evolution is consistent, i.e. to the same order, if n-loop 
beta functions are used together with matching conditions evaluated at the {n — l)-loop 
level. Firstly, we run Q;s(/i) at the one-loop order, eq. (0), with matching conditions at tree 
level, i.e. taking ai^^th) = oi^il^th) with Hth around m^. After that, we calculate as{Mz) 
by running as(/i) with two- loop beta functions and imposing matching conditions at the 
one-loop order, eq. (|l5), but taking only the first two terms on its right-hand side. And 
finally, we evaluate as(M^) according to the three-loop evolution, eq. (|^), with matching 
conditions at two- loop level, eq. (p^S]). 

We show the final results in fig. 1. We can clearly see that, as expected, the variation 
of the final prediction on a^{Mz)i as we vary the matching point around the bottom quark 
mass, is of the same order of magnitude as the next-order corrections; for three-loop beta 
functions and two-loop matching conditions it is practically flat. For comparison purposes 
we also give the error bar induced from the error in as{Mr). Given the level of accuracy, 
two-loop beta functions and one-loop matching conditions seem to be good enough for all 
purposes. 

In the preceding section we directly used the matching equation to evaluate a^^i^Mr) in 
terms of 0:3 (M,-). We could proceed in that way because the mass of the c-quark and the 
mass of the r are not so different; the logarithms in the matching equation are therefore 

■^The value of a2,{Mr) has been extracted at the two-loop level from hadronic r decays in We 
took their result directly. 

Our starting point for the quark masses are the so-called Euclidean masses |Q, — 4.23 ± 
0.05 GeV and Mf = 1.26±0.02 GeV, from which we extract the WS masses rrib = mbirrib) = 4.3±0.2 GeV 
and rric = mc{mc) = 1.3 ± 0.2 GeV. 



not large. Alternatively one could try to run a^lMr) until some intermediate scale Hth 
around the charm threshold. Then, impose eq. ([ISD with rUg = rric to get a^lfith) and 
run it until the bottom-quark threshold. This time, since we are interested only in the 
error induced by crossing the charm threshold we will use eq. (|18[) with = mf, and 
fJ'th = mb fixed to obtain a^i^mf,). Finally we run 05 (/i) from until Mz- Of course 
this procedure should give, within the level of precision of the order considered, the same 
result as before. In fig. 2 we give 0:5 (M^) as a function of the matching point fxth taken 
for the charm threshold. Although now the result depends on the matching scale ^th, 
this dependence is always a next-order correction as long as the matching conditions 
are implemented correctly. Clearly this procedure is potentially very dangerous since 
an incorrect use of matching conditions could lead to a false strong dependence on the 
matching scale. A similar consideration could be applied to the bottom quark threshold. 
Then, probably the safest procedure to run a^^Mr) until the Z mass would be to use 
first eq. ([I8|) with rUg = rric and fith = Mr to get 04 (M,-) in terms of as (Mr), then use 
again eq. (p!8|) with rUg = nib and fith = Mr to get a^^Mr) in terms of a^^Mr). Finally 
we should run a;5(/x) from Mr until Mz with the full five-quark renormalization group. 
This procedure is justified since the masses of the 6-quark, c-quark and r-lepton are not 
so different as to spoil the validity of the matching equation. Working in this way we 
arrived at the valueQ 05 (M^) = 0.123 ± 0.004, in complete agreement with our previous 
result. 

To conclude, we would like to remark on the following points: 

• Only in MO-like schemes, where Appelquist-Carazzone is realized in both beta 
functions and coupling constants, the strong coupling constant as{f^) is continuous. 
In MS-\ike schemes one should build a low-energy effective field theory and write 
scale-dependent matching conditions in order to connect the parameters of the the- 
ories on both sides of the threshold fith- Then, for general values of fith the couplings 
are not continuous although in the case of only one coupling constant it is always 
possible to find a particular fith that makes the coupling continuous. 

• Evolution is consistent, i.e. to the same order, if the evolution of the gauge coupling 
constant at the n-loop order is accompanied by matching conditions at the [n — 1)- 
loop level. 

• Different choices for the trace in Dirac space, i.e. Tr{/} = 4 or Tr{/} = 2^/^, give 
rise to different trivial modifications of the MS scheme with quite different matching 
conditions. Should one insist on having a continuous coupling across thresholds, it is 
clear from the discussion that the precise matching point will depend on the choice 
for the trace in Dirac space. For instance, working with two-loop beta functions one 
should take 

^th = nib if Tr{/} = 4 

^ItH = V2mb if Tr{/} = 2^/^ . 

By running the strong coupling constant from the Mr scale to the Mz scale, we have 
checked explicitly that the final answer is not sensitive to the exact value of the matching 
point, fith, used in crossing the nib threshold as long as the right matching conditions are 
consistently taken into account (fig. 1). Similar considerations apply when crossing the 
rric threshold (fig. 2). Finally we have shown that the correct result can be obtained by 

^We include the error induced by the errors in the quark masses, which is about 0.001 in a^{Mz). 



using the matching conditions to find 0:5 (M^) in terms of a3{MT-) and then run it with 
the full five-quark renormalization group until the Mz scale. 
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Figure captions 



Figure 1: Strong coupling constant at tlie Mz scale, obtained by running the coupling 
from its value at the M^- scale {a^^MT-) = 0.36 ±0.03), as a function of the matching point 
taken to cross the 6-quark threshold. The long-dashed line is obtained by using one- 
loop beta functions and tree-level matching conditions. The dashed hne is obtained with 
two-loop beta functions and one-loop matching conditions, and the solid line is obtained 
with three-loop beta functions and two-loop matching conditions. Error bars on the final 
three-loop result are given for comparison purposes with other as{Mz) results. 

Figure 2: Same as in fig. 1, but varying the matching point around the c-quark threshold. 
The matching point for the 6-quark is now fixed at nib- 
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